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1. Modeling of the LST System with the CMG Control Loop 

The purpose of this section is to study the dynamics of the LST control 
system, and to arrive at a simplified model for computer simulation without 
loss of accuracy. A dynamic model for the LST control system with the 
complete CMG control loop is shown in Figure 1-1. Other than the CMG 
nonlinearity model, this system is essentially the same as that in reference 
[1]. It should be noted that the limits on the current command, the gimbal 
rate command, and the torque amplifier voltage, have been neglected. The 
model also assumes a second- order pure inertial vehicle dynamics. The 
cogging torque, ripple torque, and the tachometer ripple of the CMG are 
also neglected. 

Since the LST control system has a sample-and-hold in the forward path, 
discrete-data system theory should be used when studying the system's 
behavior. However, because of the complexity of the CMG nonlinearity, an 
analytical study of the entire system is quite complex. The present task 
calls for two approaches to the analytical study of the stability of the 
system. The first analysis involves the assumption that the sampling rate 
of the system is high, so that the system may be approximated by a continuous 
data system. This is accomplished simply by shorting out the sample-and- 
hold in Figure 1-1. The second approach involves the insertion of a 
fictitious sample-and-hold at the input of the nonlinearity. This way, 
the system has two synchronized samplers, and discrete- data techniques such 
as the discrete- describing function method may be used. 

Characteristic Equation and Eigenvalues of the Continuous- Data L5T 

In this section the LST system modeled in Figure 1-1 is considered to 
be of continuous-data, that is, without the sample-and-hold. If we 
represent the input-output relation of the CMG nonlinearity between B a nd 
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T G p by an equivalent gain N, we can derive the characteristic equation 
and determine the eigenvalues of the linear portion of the system. 

The characteristic equation and the eigenvalue studies of the system allow 
us to gain understanding of the dynamic behavior of the system, and 
eventually lead to the simplification of the system model. 

With reference to Figure 1-1 # the characteristic equation of the 
continuous model of the system for stability analysis may be obtained 
from the following equation: 

1 + E Loop gains =0 (1-1) 

Or, 

1 + NO + G 4 G 5 )|§- + H 2 G 6 G s + 

+ GjG^GgGg + G-jGgGjjG^GgG^ + GgGgH^Gg = 0 (1-2) 

Dividing both sides of the last equation by the terms which do not contain 
N, we arrive at the equation which is of the form, 


where 


1 + NG(s) = 0 


G(s) = 


_ G 6 (1+G 4 G 5^ 


SA, 


and 


A 0 = 1 + H G 6 G 8 + + G 4 G 5 + G 3 G 4 G 5 G 6 

+ G 1 G 2 G 3 G 4 G 5 G 6 G 7 K t + g 6 g 8 g 4 g 5 h 2 


G-j - Kq + K-jS 


G 2 * K„s+K t 
p I 


(1-3) 

0 - 4 ) 

(1-5) 


4 


K.s+K, 

G = J 1 
b 3 Kj-S 


c . W 1+ V> 

G 4 " iTT+T“s') 


1 


6 5;yrfT^T 


- 


_ 1 


J G S 


G-7 = 


_ H 


J v s ' 


(1-6) 



Z+d oa s+k oa 


following system parameters 

are given: 

H = 600 

ft- lb- sec 

CMG Angular Momentum 

V 2 - 1 

ft-lb-sec^ 

Gimbal inertia 

*0A “ 2 - 1 

ft-lb-sec 2 

CMG output axis inertia 

d oa = 2 ° 

ft-lb-sec 

CMG output axis damping 

k oa = 1q6 

ft- lb 

CMG output axis spring constant 

K q = 5.75835X10 3 

vehicle controller coefficient 

K| = 1 . 37102X1 0 3 

vehicle controller coefficient 

K T = 2.5 

ft-lb/amp 

Torque motor sensitivity 

K b = 3.4 

volts/rad/sec 

Torque motor back emf constant 

R m = 4 ' 4 

ohms 

Torque motor armature resistance 

T e = 6.4 

milliseconds 

Torque motor armature time constant 

Kp = 216 

ft T lb/rad/sec Gimbal rate loop proportional gain 
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Kj = 9700 ft-lb/rad. Gimbal rate loop integral gain 

Ken 55 1.6X10^ volts /sec/ amp Current loop Integral gain 

T-j = 6.4 milliseconds Current loop lead time constant 

T 2 = 0.16 milliseconds Current loop lag time constant 

J y = 10 5 ft- lb-sec^ Vehicle inertia 

The numerator of Eq. (1 -»2) may be regarded as the characteristic 
equation of the overall system. Let us consider that the CMG is without 
the nonlinear friction chearcteristics. The overall LST system is linear, 
and the characteristic equation Is given by 

* 0 3 0 (1-7) 

Substituting Eq. (1-6) into Eq. (1-5), and simplifying, the numerator 
polynomial of Aq when equated to zero gives the characteristic equation. 
Thus we have 

A = b^s + b-jS** + b 2 s^ + + b^s^ + b^s^ + b^s^ + b^s^ 

+ b 8 s + bg = 0 (1-8) 

where 

b 0* R m J G J V T E T 2 I 0A 
b l = R m J G'V D 0A T E T 2 + W T E + T 2^ 
b 2 * R m J G J V Cl 0A + D 0A (T E + V + K 0A T eV + K CLI J V J G T 1 I 0A 
+ WvVoa + %£ 

b 3 ~ R m J G J V*- D 0A + ^ T 2 + T E^ K 0A^ + K CLI J V‘V T 1 D 0A + *0A^ 

+ K B W T 2 D 0A + ! 0A^ + " Z W J 2 + V + J V K CLlVoA K P 

b 4 ° R m J G J V K OA + K CLI J vV D 0A + K 0A T 1 ^ + K B K T J V^ D 0A + K 0 aV 
+ + J V K CLI^ T 1 D 0A + WV + K lW + WA 
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b 5 ' ■* C CLI J V J G K OA + VrVoA + j v k cli*- k p^ d oa + k oa t ^ 

+ M T l D OA + + K I K GLI HT lWl + K CLI J V h2 

b 6 ° j v k cli*- k oa k p + K I^ D OA + K OA T l^ + K I K CLI H ^ T l D OA + ! OA^l 
+ Woa ] 

b 7 = J V K CLI K I K OA +[IC 0 (T l D OA + *OA^ + K l^ D OA + K OA T 1^ 

b 8 = k i k cli h [ k oa k i + k o( d oa 

b 9=VcLI H VoA (1-9) 

Substituting the values of the system parameters into these coefficients, 
we have, 

b Q - 1,98697 
b 1 « 12747.9 
b 2 = 0.483289X1 0 8 

t> 3 = 0.1 92517X1 O 11 

b 4 = 0.2731 87X1 O 1 4 

b g = 0.619189X1O 16 

b g = 0.44553X10 18 

b ? = 0.163408X10 20 

b g = 0.131 112X10 21 

b 9 = 0.5362 17X1 0 21 

The eigenvalues of the system, or the roots of Eq. (1-8) are found 
and are tabulated below: 
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-4.56811 + j4. 68534 ) 

V vehicle dynamics 
-4.56811 - j4. 68534 J 

-1563.11 CMG torque motor Armature time constant 

-39.8505 + j 40. 9207 
-39.8505 - j 40. 9207 
-12.4647 + j 748. 038 

V CMG output axis torsion dynamics 
-12.4647 - j748.038 J 

. -3072.85 + J3554. 86 1 

7 CMG current loop controller 
-3072.85 - j 3554. 86 J 

The contribution of the sectors of the 1ST and the CMG dynamics to these 
eigenvalues are indicated. 

It is essential to investigate the importance and weight of each of 
these eigenvalues upon the dynamics of the entire system. Since some of 
these modes have very short time constants and high-frequency oscillations, 
it is expected that a digital computer simulation of the entire exact 
model will be time consuming and costly. 

We shall show that for computer simulation and analytical purposes 
some of these fast and oscillatory modes of the system may be neglected. 

Figure 1-2 illustrates the root locus diagram of Eq. (1-8) when K Q 
and K 1 are varied proportionally from their nominal values which are listed 
earlier. In Figure 1-2 the variable parameter is indicated as K which is 
the multiplier of the nominal Kq and Ky Table 1-1 gives the location 
of the roots as a function of K. 

It is interesting to note from the root locus diagram that for the 
nominal Kq and Kj both the LST vehicle dynamics and the CMG dynamics have 
a relative damping ratio of approximately 0.707. Furthermore, the 


} 


CMG dynamics 


-3073+13555 


-3073-13555 


s- plane 


-l2.47+j?48,K*0.1 l 

\ 

K*100 


-'44+|45.l, K*0.l 


K e 0«5 


-39.65 + J40.9, K“I.O 


K-2.0 


K*I0 K*9 K«5 I K-5 9 10 


-85.85, K«KX> K*9.7 -19.5 


K*0.l .. *l el K“2 


K*I00 


-I2.47-/748 




Figure 1-2. A root locus diagram of Eq. (1-8) when K Q and K-j are varied proportionally 



/K"O.I 

K« < 1.5 




Table 1-1. Eigenvalues of the complete LST as 


Kq - KKq (nominal) and K-j ~ KK^(nominal) vary. 
Kq (nominal) = 5.75835X10^ 

K, (nominal) = 1.3710ZX10 3 


r 


K 


Roots 


0.1 

-0.379 

± 3 1.837 

-1563.1 

0.5 

-2.041 

± j3.853 

-1563.1 

1.0 

-4.568 ± j4.685 

-1563.1 

2.0 

-6.295 

-19.522 

-1563.1 

5.0 

-4-680 

-64.314 

-1563.1 

9.0 

-4.441 

; 

-82.631 

-1563.1 

9.7 

-4.422 

-84.921 

-1563.1 

10 

-4.415 

-86.866 

1 -1563.1 

100 

-4.220 

-156.13 

-1563.1 


-44.034 ± j 45. 091 
-42.374 ± j 43. 335 
-39.850 ± j 40. 921 
-31.515 ± j 35. 779 
-9.943 ±j47.136 
-0.925 ± j58. 522 
0.206 ± j60.068 
0.667 ± 3*66.774 
45.437 ± J 1 29.14 


-3072.8 ± j 3554. 8 


-12.470 ± J748.04 
-12,467 ± J748.04 
-12.465 ±j748.04 
-12.459 ± j 748. 04 
-12.443 ± J648.04 
-12.422 ± J748.04 
-12.418 ± J748.04 
-12.417 ± j748.04 
-11.94 ± . J748.03 


-3072.8 ± j3554.8 
-3072.8 ± 3*3554.8 
-3072.8 ± j3554.8 
-3072.8 ± j 3554. 8 1 
-3072.8 ± j 3554. 8 
-3072.8 ± j 3554. 8 
-3072.8 ± 33554.8 
-3072.8 ± J3554.8 
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eigenvalues at -1563.1, -3072.8 ± j3554.8, and -12.465 ± j748.G4 are not 
sensitive to the change of the loop gain at all. 

Another way of investigating the significance of each of the 
eigenvalues toward the dynamics of the entire system is to evaluate the 

residue of the closed- loop transfer function that corresponds to each of 

/ 

the eigenvalues. 

From Figure 1-1 the closed-loop transfer function of the LST system 
without the CMG nonlinearity is written 

9 W G-* GrtG-GiiGrGx- G-*K T 

_J/ = 1 2 3 4 5 6 7 T (1-10) 

X A o 

After simplification, Eq. (1-10) becomes 

6 V _ l < T K CL i K i H ( a o s4+a l s3+a 2 s2+a 3 S+a 4^ , . 

X A 

where A is given, by Eq. (1-8), and 
a 0 = ^0A K 1 T 1 

a l = D 0A K 1 T 1 + W K 1 + K oV 

a 2 " >0A T 1 K 1 + °0A (K 1 + W + Wo 

a 3 = K 0°0A + W K 1 + W 

, ' a 4 K 0 K OA ' 

Substitution of the system parameters into Eq. (1-11), and performing 
partial faction expansion, the residues of Eq. (1-11 } at its poles are 
tabulated in Table 1-2. It is interesting to note that only the residues 

i 

which correspond to the CMG and the vehicle dynamics are of significance. 
The contributions from the" eigenvalues at -1563.1, -3072.8 ± j3554.8, and 
-12.465 ± j 748. 04 are very insignificant. 
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Table 1-2. 

Residues of Eq. (1-11) at 

Its poles 


Roots 

Residues 

Vehicle dynamics 

-4.56811 ± j4. 68534 

14.3102 Tj0.5256 

CMG torque motor 

armature time constant -1563.11 

0.002778 

CMG dynamics 

-39.8505 ± j40. 9207 

-14.2967 ± j 12 . 7188 

CMG output axis 
torsion dynamics 

-12.4647 ± j 748. 038 

-0.013248 ± jO. 0001275 

CMG current loop 
controller 

-3072.85 ± j 3 554. 86 

-0.001565 T jO. 002523 
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Still another method of investigating the effects of each of the 
system components on the system dynamics is to make the frequency-domain 
plots for G(s) of Eq. (1-4). These plots are also useful for the 
continuous-data describing function study of the system. In Figure 1-3, 
curve 1 gives the G(s) plot in decibels versus phase for the 9th-order 
LST system. In the same figure, curve 2 gives, the G{s) plot when the 
output axis torsional dynamics of the CMG are neglected; that is, I QA - 0, 

°0A * °> K 0A ‘ °' 


In this case, G(s) is simplified to 


G(s) - 


R m s ( 1 +T 2 s ) ( 1 +T e s ) + K CL i ( 1 +T| s ) 


(M2) 


where 


A - W2V 5 * Wg< T 2 * V s ■ 

+ ( R m'V + K CLI J G T l + K B K T T 2^ s3 + i K CLI J G + Vt 

+ k cli k p t i )s2 + k cli (K; p + k i t i )s + k cli k i 0 - 13 ) 

In addition to neglecting the CMG torsional dynamics a further 
simplification of the LST system can be conducted by neglecting the 
time constants Tj, and T^, as these are small when compared with the 
time constants of the vehicle and the CMG. Figure 1-4 gives the block 
diagram of the Gimbal rate control loop and the CMG with these simplifications 
Using the block diagram reduction technique, the system in Figure 1-4 is 
reduced to the block diagram of Figure 1-5. Since K^/R^ is very large, 
the transfer function of the inner loop involving K CLI and R m is approximately 
unity. Furthermore, since Kg/K^j is very small, the back emf loop of 
Figure 1-5 may also be neglected. The final simplified gimbal rate control 
loop and CMG dynamics are represented by “the block diagram of Figure 1-6. 


The transfer function of the simplified CMG control loop is 
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Y 


KpS+K r 

J G s 2 +K p s+Kj 


0-14) 


Substituting the system parameters into the last equation, we have 


V_ 102. 857 (s+44. 9074) 

Y ' ( s+51 . 428+J44. 43 ) (s+5l . 428-j 44. 43) 


Now replacing the CMG control loop in Figure 1-1 by the simplified 
diagram of Figure 1-6, we have the block diagram of Figure 1-7 for the 
simplified LST system without sampling. 

The transfer function for G$s) in Eq. (1-4) is simplified to 



0 g J v s 4 +0 v K s 3 +O v K i s 2 +KjHK 1 s+K i K 0 H 

(1-16) 

i Figure 1-8 shows the plot for the G(s) in Eq. (1-16). It is 
interesting to note that the three G(s) plots in Figures 1-3 and 1-8 
are very similar for frequencies below 1000 rad/sec. This proves 
that the LST system may justifiably approximated by the simplified system 
of Figure 1-7. 















21 


PRECEDING PAG® BLANK NOT FILMED 


Then can be generated by intebratlng both sides of Eq. (2-1) with 
respect to t. It Is observed that % Is easily found as a function 
of x. However, the simulation model actually has Tq F generated by a 
function generator with as the input. The simulation block diagram 
for the CMG friction nonlinearity is shown in Figure 2-1. Notice that 
this model has 0^ as input and T^p as the output. 

Referring to Figure T-1, _ the block diagram of Figure 2-1 may 


replace the blocks for the simplified nonlinearity between e G and Tgp, 
It has been demonstrated experimentally [3] that for the solid 
rolling friction the relation between T^p and T^p may be approximated 
by a square-law expression. 


dT 

T 1 s 

GF da 


y ^ t gf ■ t gfo^ 


t gf - t gfo 


(2-2) 


where Y is a positive constant. 


T GF ■ 0 t gf * t gfo 


(2-3) 


However, for the CMG friction, the frictional torque is also velocity 
dependent, as shown in Figure 2-1. Therefore, Eq. (2-2) should be written 


t gf * y ^ t gfi " t gfo^ Z ( 2 “ 4 ) 

T GFI " ( 2 ~ 5 ) 

For simulation purpose, the block diagram of Figure 2-1 can be 

easily programmed on the digital computer, and the program can be used 
for any input. 

Our objective Is to investigate the behavior of the 
nonlinearity under a sinusoidal exci tati ohT so that the describi ng 
function can be derived. In the present situation, however. 
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we shall show that the input-output relation of the CMG 
nonlinearity is different for sine and cosine inputs. It should be 
noted that the describing function is a pseudo-sinusoidal steady- 
state analysis so that either sinusoidal or cosinusoidal inputs may 
be assumed. These are shown by the following derivations. 

(1) Let ' = A sinwt (2-6) 

then 0 r = - — cosoit (2-7) 

G a) 

For 0g >_ 0, Eqs. (2-1M2-4) and (2-5) lead to 

dT fs r p 

IT = ^Asinut ^ t gf - W {2 ' 8) 

2kTt<o>t < (2k + 1 )ir, k = 0, 1 , 2, .... 

Solution of the nonlinear differential equation of Eq. (2-8) yields 

-f — Tjr*- — = - coSwt + C (2-9) 

! gf~*gfo w 


The constant C is found by specifying the initial value of T^ at t - 0; 

t = 0 T GF ~ “ T GFi T GFi > 0 ( 2 “ 10 5 

Thus, 


c = lA + ! 

* T GFi +T GF0 


(2-n) 


Substitution of the last equation into Eq. (2-9) and simplifying, we have 
the solution 


T GF = 2^“ cos ^)" R+l 
T GF0 -KTT+ irO'COSwt) 


(2-12) 


R+l 2 

which i s valid for 6 q > 0 , or 2k7r < cat < (2k + 1 )tt, k = 0, 1 , 2, . . . 
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where 


and 


b * 


R - 


2 T GF0 YA 


U) 


T 

GFi 

t gfo 


For 6g£0, Eqs. (2-1), (2-4), and (2-5) lead to 


(2-13) 


(2-14) 


dT r - 9 

“dtT = yA sinwt (’ T gp “ T qpq) 

(2k+l)TT £ wt <_ (2k+2)TT, k = 0, 1, 2* .... 
The solution of Eq. (2-15) is 


1 


' t gf' t gfo 


iA 

03 


cosoit + C 


(2-15) 


(2-16) 


The boundary condition is T GF = -T GFi at cot = 2 it and its integral 
multiples. Thus, 

c = ^ + 


“ T GFi" T GF0 


(2-17) 


Substitution of C into Eq. (2-16) and simplifying, the final solution is 


fgp ” R-l ~ 2 ^~^OS0)t) 

T GF0 j ( l-COSwt)+|j~y 


(2-18) 


which is valid for 0^ £ 0, or (2k+l )tt < ost < (2k+2)iT, k - 0, 1 , 2, . . . , 
and 0 q is given by Eq. (2-6). 

“ T Gpi » which represents the value of T^p at wt = 2ku , k = 0, 1, 2,..., 
can be determined by matching the values of ^p/Tqpq at cot = (2k+l)7r, 
k ~ 0, 1, 2, ...» for the two ranges of 0^. Setting cot - tt and 
equating Eq. (2-12) to Eq. (2-18), we have 
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b - 


R+l 


R-l 


- b 


1 

Rfl 


+ b 


b + 


(2-19) 


R-1 


which is simplified to 


R Z + gR-1 = 0 


Solving for R gives 
R 


■ • 


?+T 


(2-20) 


( 2 - 21 ) 


Since T^ £ T G p G and Tgp. > 0, R _< 1; the plus sign must be chosen in 

Eq. (2-21). 


R 


V i | M 1 ■' — » 

b 2 +l _ T GFi 

~ 3 ~~ T^r 


( 2 - 22 ) 

b‘ 1 GFO 

The same result for R can be obtained by use of the boundary condition at 
wt - 7r which is Tgp = T^p. 

The consequence of defining e Q as a sine function, Eq. (2-6), is that 
Tgpj 1S dependent on w, since b depends on to. Figure 2-2 gives’ the 
T GF^ T GF0 versus 0 g^ A relations for various values of A/w. Notice that 
the maximum value of Tqjt/Tqpq is different for different ratios of A/ou. 
Figure 2-3 illustrates the relations between T GF /T GFQ and e^/A with 
Qq 9iven by Eq. (2-7). Note that the e G axis is normalized by W /A, so 
that as w varies not only the bounds of T Gf /T GP q will vary but also the 
bounds of Oq. These curves are also plotted for various values of A/to 
indicated. It would seem that by defining e G as a sine function as in 
Eq. (2-6), there is an advantage of using the nonlinear characteristics 
of Fig. 2-2 for analytical purposes. In principle, either e r or e„, 
but not both, may be used as the input to the nonlinearity , provided that 
the corresponding nonlinear characteristic is used. 


GFO 



Figure 2-2, Normalized frictional torque versus Qq/A for CMG nonlinearity, with sine function input. 


ro 

o> 





- 0.8 - 0.6 - 0.4 - 0.2 0.0 0.2 0.4 06 - 0.8 1.0 


0 g u>/A 

Figure 2-3. Normalized frictional torque versus e^/A for 
CMG nonlinearity with sine function input. 
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This is because 9 q and 0^ are directly related'. In other words, the 

T GF^ T GF0 versus 0 g characteristics contain the same information as 
those of T G p/TgpQ versus 6g. 

(2) Let 0g - A cos wt (2-23) 

Then 8^ = -Awsin wt (2-24) 

For 0 G < 0,-Eqs. (2-1), (2-4) and (2-5) lead to 

dT rF ~ 

dt *\ yAtosincot (“Tgp “ ^qF(P (2-25) 

2k7T < (2k+l)7T, k = 0* 1 , 2 , . . . 

The solution of Eq. (2-25) is 


-1 

T +T 
GF GFO 


- yA cos wt + C 


(2-26) 


At t = 0, the boundary condition on T 
Equation (2-26) gives 


1S T GF * T GFi ’ where T GFi > ° - 


C = - yA - 


T GF + T GF0 


(2-27) 


Substituting Eq. (2-27) into Eq. (2-26) and simplifying, the complete 
solution is 


Tgp {^4.'] ~ 2"(l-cosa)t) 

T GF0 |{l-coswt)+Jy 

which is valid for 2k?r ± ufc < (2k+l)?r , k = 0, 1 

a * 

For 0 q > 0, 

"dt = " YAa)sinwt(Tgp - T^) 2 
(2k+l )ir £ o)t < (2k+2)ir, k - 0, 1 » 2, ... 


(2-28) 

2, — , or 6 q £ 0, and 
(2-29) 


(2-30) 
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The solution of Eq. (2-30) is 


-1 


T GF ~ T GFQ 


= yAcoswt + C 


(2-31) 


The boundary condition is T GF = T GFi at ait * 2tt. Thus, Eq. (2-31) gives 


C = -yA - 


T GFr T GF0 


and the complete solution is 


GF R-l + 2* 

a 1 

2^1-COSO)t) + g-y 


GFO 


(2-32) 


(2-33) 


which is valid for 'e Q > 0, or (2k+l)7r < tot < (2k+2h a k » 0, 1, 2 
Matching Eqs. (2-28) with Eq. (2-33) at wt = 7 r, we have 
2 2 

R + f- R - 1 = 0 
a 


whose solution of R is 


R * -i + 



, 2 T 

a +1 _ 'GFi 

~ T~ " ? — 
a *GF0 


(2-34) 


The significance of this result is that T rc ./T rr . 

brl bru 


(2-35) 

is not a function of w. 


since a depends only on y» A, and T^q. It is interesting to note that 

Eq. (2-22) and Eq. (2-35) have the same form, except that in the former 

equation b is a function of w, whereas in the latter equation a is not. 

Figure 2-4 shows the ^p/T^g versus Qg/A characteristics for 

several values of A. The maximum value that T GF reaches for all values 

of oj is the same. Figure 2-5 gives the T GF /T GF0 versus 0 Q /A 

characteristics with the © G axis normalized by Aw. From these curves, 

it would seem that there is a clear cut advantage in using the T 0r /T 0 

GF GFO 

versus e Q /A characteristics with 0 Q = Acoswt as the input. Referring 
to Figure 1-1 , for describing function analysis, one difference between 



Figure 2-4. Normalized frictional torque versus e G u/A for 
CMS nonlinearity with cosine function input. 


CJ 

o 




0 G a Aco$ (wt) 

8 q =*-Aaisin(wt) 


Figure 2-5. Normalized frictional torque versus e G /A for. 

CMG nonlinearity with cosine function input. 




the use of 8^ and 0^ as the input to the nonlinearity is the linear 
transfer function G(s) which the nonlinear element sees. 

Table 2-1 gives a comparison of the expressions of 
for the two types of inputs. 'It is interesting to note that if A/to 
in b is replaced by A, the expressions for T G p/T G p 0 differ only by a 
sign for the two inputs. A similar effect is achieved by replacing 
-A/w by A in Eq. (2-7) for e Q , and A by -At* in Eq. (2-6) for e G . 

With the above changes , , the torque expressions for the two types of 
inputs become identical. This means that the curves in Figs. 2-2 and 
2-5 and those of Figs. 2-3 and 2-4 are identical for the corresponding 
values of A/w and A, respectively. 
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A More General CMG Frictional Nonlinearity Model 

A more general expression which is valid for any input can be 
derived for the CMG frictional torque characteristic* We start with 
the square-law expression in Eq. (2-4), 

dTpp p 

T GF = de G " y ^ T GF1" T GF0^ T GF - T GF0 (2-36) 

where > 

T GF1 = ^GfSG^( 0 g) (2-37) 

Equation (2-36) can be integrated with respect to e G directly. 

The results are: 


T = ^ 9 G~ 9 Gi )r(T GFi' T GF0^ T GF0' T GFi 
GF ^ S G” e G1 ^ T GF1 ' T GF(P ’ 1 


e G — 0 


: 9 G~ 9 Gi ^^ T 6F’i + T GF0 ^GF0~ T GFi 
( e~- )Y(T« 1 -.+Y nf . rt ) - T 


G Gi 7 Y v GFi GFO 


9 G^° 


where 


= initial value of 
T GFi ~ value of 

These expressions for the frictional torque are. valid for any arbitrary 

input e G and initial conditions at the beginning of the process. 

* 

When changes sign during the process, the proper expression should 

be used, and the initial conditions should be appropriately matched at 

the switching point. It is simple to show that when e« * - ^cosat, 

b to 

e Gi = ~ u’ Eqs ‘ (2 * 38 ) and ( 2 ‘ 39 ) revert t0 Eqs. (2-12) and (2-18), 

■ respectively. Similarly, if e G = Acoswt and e.. = A, Eqs. (2-38) 
and (2-39) become identical to Eqs. (2-28) and (2-33).. respectively. 
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3- A Continuous Describing Function for the Gimbal Friction Nonlinear ity 

In chapter 2 the nonlinear characteristic of the gimbal friction 
is analyzed with sinusoidal inputs. It seems that there is an advantage 
of using a cosine wave as the Input, since, then, the fri cti on a! torque 
saturation level is not a function of frequency. For analytical 
purposes, the Tgp/Tgpg versus e^/A curves of Figure 2-4 seem to be the 
most convenient to use. 

In the last bimonthly report. No. 1-73, the gimbal frictional 
nonlinearity is approximated by a four-sided polygon. In this section 
a better approximation is made by the six-sided polygon shown in 
Figures 3-1 and 3-2. It is apparent from Figure 2-4 that when the 
input amplitude changes the slope k of the nonlinearity also changes. 

The continuous-data describing function for the nonlinearity in 
Figure 3-1 and 3-2 is now derived. 

The input to the nonlinearity is a cosine waveform shown in Figure 
3-1 and given as 

0g - Ecoswt (3-1) 

The output waveform of the nonlinearity as a function of .cut is shown 
in Figure 3-1. The four parts of this waveform are described by the 
following set of equations. 

f(e 6 ) = kEcos# - kE 0 < $ <_ ai 

< f 

s k£cos<ji + k£ 7T <_ <*2 

= M - a 2 < <{> <. 2tt (3-2) 

4> - a)t 

a-j =' cos^O-M/kE) 

a 2 33 a>j + tt 


where 


(3-3) 

(3-4) 

(3-5) 
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The Fourier series for the waveform of the output of the nonlinearity 
is given by 

f(t} = Aq/ 2 + A^coswt + BjSinwt 

+A2cos2cot + B2Sin2o)t + (3-6) 

where the coefficients are given by 


A^j = 1 /tt j ■. f(t)cos(Nwt)dwt ( 3 - 7 ) 

j2tt 

B n 35 1 /tt j f(t)sin(Nwt)da)t (3-8) 

N *0, 1,2, ... 

Since there is no bias on the output signal from the nonlinearity, A Q - 0. 
If the higher frequency harmonics, N « 2, 3, are neglected, the 
output waveform is given by 

f(t) -J^f cos(<j>- tan" 1 B 1 /Aj) (3-9) 

The phasor notation of f(t) is 


F(jw) A* + expC-jtan^B^A^ (3-10) 

The describing function of the nonlinearity is defined as 

N(E) = F(jw)/E(ju) /E-expC-jtan^Bj/Aj] (3-11) 

= A ,/E - jB-j/E 

By using Eqs. (3-2) through (3-5) the coefficients in Eqs. (3-7) 
and (3-8) are solved for N = 1 to yield 

Aj = kEcos'^O - M/kE)/7r 

+ J2 MkE - M 2 (M-kE)/frkE (3-12) 

B 1 = -4M/w + M 2 /^ (3-13) 
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Describing Function Plots 

When evaluating the describing function it is convenient to define a 

\ L . 

normalized input amplitude H as 

H - 2kE/M (3-14) 

Using equation (3-14) in Eqs. (3-12) and (3-13) and dividing both sides 
of these equations by E gives 

A-j/E '» kcos _1 (l -2/H)/« +4k(T - H/2)VH^/H 2 ir ■ (3-15) 

B^E = -8k/Hir +4k/H 2 ir (3-16) 

The above equations show that the describing function will depend 
only on H and k. 

Figure 3-3 shows the magnitude versus phase curves for -l./N(E) 
for different values of the slope of the' nonlinearity , k. In these 
curves the magnitude of the normalized input signal to the nonlinearity* 
H, is varied from 1.0 to 100. When the magnitude of H becomes less 
than 1, the gimbal torque characteristics decrease as shown in Figure 
3-2. For these amplitudes of H, M is set to M * 2kE and the magnitude 
and phase of l./N(E) remains constant as H is varied. This can be 
seen by substituting M =* 2kE into equations (3rl2) and (3-13). 

Aj/E = k 

Bj/E - -4k/ir 


(3-17) 

(3-18) 
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Figure 3-3. Describing function plots of CMG frictional nonlinearity. 
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4 * Computer Simulation of the Simplified LST System 

Since the describing function analysis of the CMG nonlinearity has 
been given in the last chapter, a computer simulation of the LST system 
is necessary to corroborate the analytical results. It should be noted 
that the describing function analysis is carried out only for the 
straight-line approximation of the nonlinearity as illustrated in 
Figure 3-1. The simulation model of the nonlinearity also has the same 
characteristics. Future work will include the derivation of the 
describing function of the exact nonlinearity as modelled in Figure 2-1, 
as well as a computer simulation with the exact model. 

The simplified LST system is represented by the block diagram of 
Figure 1-7. The linear transfer function which the nonlinear element 
N sees is given by Eq, (1-16). The frequency- domain plot of G(s) is 
plotted in Figure 1-8 in db versus phase coordinates. The portion of G(s) for 
a) = 2.5 to 6 is given in Figure 3-3. For k * 50,000 ft-lb/rad, the 
G(s) curve intersects the -1/N(E) curve at a frequency of 5.1 rad/sec 
and H - 9. The saturation level of the gimbal frictional torque is set 
at M = 0.078 ft— lb. Therefore, the corresponding amplitude of oscillation 
of the gimbal position^, is 

E « MH/2k * 7.2 x 10" 6 rad. 

For the computer simulation, the input to the LST system, x » is 
set to zero, along with all the initial states, except for the vehicle 
position e v . The initial value of e y is set at 5 x 10‘ 5 rad, which is 
chosen so that the input signal to the nonlinearity, e^, would be large 
enough to cause the torque to saturate, while at the same time the 
limiting value of the input signal given in [1] is not exceeded. 
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To fully understand the results of the computer simulation of the 
1ST system it is useful to represent the computer modeling of the 
nonlinearity by the flow diagram of Figure 4-1. Due to the inherent 
memory of the nonlinearity it is necessary to initially define the 
position of the last switching point and the last value of the velocity 
before entering the flow diagram. 

The following quantities are plotted from the simulation runs: 

6 y s vehicle position (rad) 
u> v * vehicle velocity (rad/sec) 

0£ = Gimbal position (rad) 
o) G = Gimbal velocity (rad/sec) 

Tgp ® Torque output of the nonlinearity (ft-lb) 

Error = Error input command (rad/sec) to the CMG 

= x " K 0 6 v K l w v 

i 

Figures 4-2 through 4-7 show the plots of the above listed quantities. 
It may be noted from the plot of T 6F in Figure 4-6 that the system has 
a sustained oscillation. This oscillation is not seen on the other plots 
because of the large initial transients., Figures 4-8 through 4-13 show 
the continuation of Figures 4-2 through 4-7 with proper scales. These 
figures show that the amplitude and frequency of the sustained oscillation 
are quite close to the values predicted In Figure 3-3. 

The simulation results show that the LST system is stable for smaller 
values of k, such as k * 10,000. 


ER 
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RETURN 


Figure 4.1. Computer flow diagram of CMG nonlinearity. 
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5. A Describing Function of the CMG Nonlinearity 
Using the Analytical Torque Equation 

A describing function of the CMG frictional nonlinearity was 
derived earlier using the straight-line approximated input-output 
relation between the frictional torque T^p and the CMG angular dis- 
placement 

However, it is possible to derive a describing function for the 
CMG frictional torque using the analytical relation between T^p and 
It has been established that the frictional nonlinearity of the CMG 
can be described by the square- law relation. 


dT GF _ ,, , ,2 

dSg GFI " GFCr 


(5-1) 


where 


T GFI * t gf sgn <V 


(5-2) 


Tqfq a saturation level of T^p 
Y = positive constant 

Carrying out the integration on both sides of Eq. (5-1) yields 


0 G + C 1 ’ 


-1 


y ^ t gf ' t gfo^ 


S G — 0 


(5-3) 


e G + C 2 


-1 


y(t gf + t gfo j 


0 G i° 


(5-4) 


where and C 2 are constants of integration, and 
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T GF e G > 0 


(5-5) 

t gf ®g i 0 


(5-6) 

given by 



1 

Gi y(T + - T 

Y ' GFi 1 GFO) 

®G - 0 

(5-7) 

1 

Y ^ T GFi + T GF0^ 

®G 1 0 

(5-8) 


where 6 q. and Tgp. denote the initial values of e G and Tgp respectively, 
For a sinusoidal input, 6g is represented by 

0 G = A cos cot ( 5 _ 9 ) 


It is important to note that for the input of Eq. (5-9) 0 Gi = -A when 
0 q ~ and e Q-j = A when 0g <_ 0. 

Solving for ftp and T~ from Eqs. (5-3) and (5-4), respectively, \ 


GF y(A cos wt + C- 


e G i 0 


(5-10) 


9 G £0 


(5-11) 


C| = A- 


Y ^ T GFi " T GF0^ 


(5-12) 


C 2 = -A - 


Y ^ T GF1 + T GF0^ 


whe re 
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a = 2yAT GF0 


(5-14) 

(5-15) 

(5-16) 


With the describing function method, the frictional torque T 6fr may be 
approximated by the fundamental component of the Fourier series. The dc 
component is zero, since the input-output relation is symmetrical 
about the zero- torque axis. 

Thus, 



sin a)t + B.] cos cot 


Substitution 


\ + ^1 ' cos 


(p = tan 


-1 A 1 


i r 2lT 

viL 1 


T rtr sin at diot - - 
GF 7T 


1 Pv 


GF 


i fy 

* K GF 


2ir 


B, = T rr COSojt dwt = . . , n _ 

TT ( 0 GF TT J 0 GF 


- 1 r 
" T Jn 


' + 

TT 


i f 2 y 

/ gf 


of Eqs. (5-10) and (5-11) into Eq. 


sinwt da)t 


sinwt dtot 


cos tot dwt 


cos wt dtot 


(5-19) gives 


(5-17) 

(5-18) 


(5-19) 


(5-20) 
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_ 1 


-1 


A 1 = ? ] n YTESit + Cj s1na * du)t - ? j 




T G p 0 siriojt dwt 


r2ir 


+ ¥ L vTF^LrTTT) s1n “ t dut + ¥ j T GF0 sinut dwt (5 ~ 2 ^ 

" ' I *-7r . 


Thus 


A 1 - ^ £n (A cos ait + C 2 ) 


+ ttAy cosa)t + c ]) 


7 r 2T 


GFO 


7T 


0 

2tt 2T 


GFO 


u 


TT 


A 1 ~ ttAy 




/C 2" A \ 

v ^nt) + ^ 


1 C-j+A) 

■ c r A ^ 


4T 


GFO 


(5-22) 


(5-23) 


In arriving at the last expression It Is noted that 
C 2 < -A 

C 1 > A 

and C-j = -C 2 over their respective ranges of ® G> Equation (5-23) is 
simplified further to 


4T 

A-j « - 
l TT 




|'(C 2 -A)(C 1 +A) 

r^+AlTc^I)] 


4 7 S ^i +A 


it GFO tiAy [C r Aj 
Now substitution of Eqs. (5-10) and (5-11) into Eq. (5-20) yields 


(5-24) 


I- 1 / 


i 

q yTa coswt” C 2 ] 


COStvlt d(i) t 


ir 

•■L 


TgP0 cos-wt doit 


if. 


2 71 __-l i f 2lT 

y[Acbsa)t + C, ] C0Swt ^oit + ^ f t qjtq coswt dwt (5-25) 
1 •'tt 
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Evaluating each of the integrals in the last equation, we have 


i 


rZ'H 

TqFO C0SCjjt d(A)t * Tgco cosoit dwt = 0 
0 J 7T ” u 


(5-26) 


? p 

Since < -A, C 2 is always negative, and > A , the first integral of 


becomes 


or 


_ i r -i i 

1 1 ^ tt L y{Acoswt+C ? ) coswt dwt = ^ 


-tot 

A I 


7T 


0 


A f 

WA j n 


dtot 


Acoswt+C 




1 2 

l 1 ~ " YA + titA ^ 


tan 


,](C 2 -A)tan((pt/2) j 


SE7 "" 


TT 


0 


(5-27) 


yA yA^-A 2 


(5-28) 


Where the fact that C 2 is negative has been used. Also, tan ir/2 is 
taken to be +•» since ut/2 expands from 0 to tt/2. 

Similarly, the third integral of B ] in Eq. (5-25) is written 


= if 

* TT 


l 2 =^ 


-1 


Y(Acosaft+C-| ) C0Suit dut 


J.+ Af 2 

YA TtyA 




tan 


(C 1 -A)tan( &J t/2) ] 

vip 


2tt 


'1 


4 + 

yA \T? F 
yA V Cj -A^ 


(5-29) 


In arriving at the last equation, we have recognized that Cj > A and 
have used that tan tt/ 2 = since in this case to t/2 expands from 77/2 to tt. 


Thus , 
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B 1 - h + ! 2 


2 C 2 + C 1 

Y ya/c|-a 2 Y A^jc 2 -A 2 


__2 

yA 





The describing function in complex form is written as 


(5-30) 


N(A) = (5-31) 

where A-| and B ] are given by Eqs. (5-24) and (5-30), respectively. A 
digital computer program for the computation of N(A) and -1/N(A) is listed 
in Table 5-1. The constant A is represented by E in this program. The 
parameters of the nonlinearity are: 

Tqpq “0.1 ft-1 b 
Y = 1.38 x 10 5 


Figure 5-1 shows the magnitude (db) versus phase (degrees) plots of 

-1/N(A) for y = T. 38 x 10 4 , 0,69 x 10 5 , 1.38 x 10 5 , 0.69 x 10 6 , and 
6 

1,38 x 10 , as the magnitude of A varies. Note that as A becomes 
large, the magnitude of -1/N(A) approaches infinity and the phase 
approaches -270 degrees. As A decreases, the magnitude of -1/N(A) decreases 

i , 

and the phase approaches -180 degrees. In the limit as A-K), -1/N(A) 

2 

goes to “1/YT GF0 * 


Asymptotic Behavior of -1/N(A) for Very Small Values of A 
Figure 5-1 shows that as A approaches zero, the magnitude of 
— 1 /N ( A ) in db approaches 20 1 o9 10 [1/yT gf() 2 ] and the phase is -180 
degrees. The asymptotic behavior of -1/N(A) for very small values 
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Table 5-7 


L 



.100 " LST CONT IlidUS DESCRIBING FUNCTION - EXACT CM6 HONL I NEAR I TV 

1.000 COMPLEX GV>6N^16 

1 .500 REALMS P < £ 0 > « P I * RAD * TO * GAMMA *E ST ART i» E V A AS R • TGF I y T G FN y T6F P * C 1 ,C2 
!*A1 >B1 

1 .600 REALMS AZ ‘4 ■ 

£.000 P 1=3. 14159 >4. 

3.000 RAD=1 80 . - PI 

4.000 to= . i ' : 

4.500 s=io. 

5.000 GAMMA=SM .38E5 ■* V 

6.000 ESTAPT=1 .E-13 ' 

7.000 NP=£ -'.y ; 

8.000 ND= 15 ' ^ 4 

9.000 P(l) = l . ' ^ " 

10.000 P<£>=5. ... 

10. £00 WRITE <6 » 1 00) 4 4: 

1 0.400 WRITE <6 *101) , V \J ' 

1 1 .000 BO 1 J=i »NB 4’'. 

12.000 DO 1 I = 1jNP ' V *’- v . 

13.000 E~E$ThRT^P< I }♦< 1 0 . ♦♦< J-l ) > , . 

14.00 0 AR=£ .♦GAMMR^E^TO . 

15.000 R~<-1 . -AB ) +BSQRT < < AA*BN+T . > < CAB ♦AN) ): 

16.000 TGFI=R^TO 


17.000 TGFN=TGFI 

18.000 TGFP=-TGFI 

19.000 C1=E-1 < GAMMAS CTGFP- TO > > \ 

£ 0 . 0 0 0 C6--E-1 . / < GAMMAS < TGFN+TO > ) ’■ 


£1.000 A 1 = < -4 . ♦TO/P r > ♦ < 1 . ✓ C P I ♦GRMMA^E > > ♦DkQG < < < C 1 +E > ♦ < C£-E > > v < c. C 1 -E > ♦>. 
£1 .100 CC£+E>)> J 

£1.500 RZ=DLDG< < <C1 +E)^CC£--E) )/< <C1-E>^<C£+E) ) ) 

£1 .600 A 1 * < -4 . ♦TO^'P I > + CRZX CPI ♦GAMMA fE > > 

££.000 B 1 = < - 1 ,/< GAMMA^E ) ) ♦< £ .♦e£^BSQRT<C£4C£~E^E)-“Cl^BSQRT<Cl^Cl-E^E> > 
££.100 A1=H1,-E ■ * • 

££.£00 B1=B1 'E • 

£3.000 GN=BCMPkX<Bl »-Al> 

24.000 GV=-1 .>GN . ' 

£5.000 G1=REAL<GV) 

£6.000 6£=AIMAG(6V) ' ' 

£7.000 GMAG=CABS<6V> 

£8 . 000 GBB-£0 .♦AkOGl QCGMAG# 

£9.000 GPHASE=RAD«ATAN2<G£>G1) 

3 0 . 0 0 0 I F C GPHASE . 6E . 0 . > GPHASE«GpHASE-36 0 . 

33.000 WRI TEC6 >1 0£>E>TGFI rGPHBSE >G1)B »GMBG 

34.000 1 CONTINUE 

35.000 100 FORMATS CONTINOUS DESCRIBING FUNCTION FOR CMG NONL I NEAR I 

TV •- > . . ■% ■' ' •. • 

36.000 101 FORMAT (V *8X s ''E y j 1 1 X » 'TGFI v > 1 OX PHASE 'Vs 1 OX > ' BE ' >9X ? 'MAGNTTU 

BE'V> ' ■■ ■ o • 

37.000 1 0£ FORMAT ( 1P5E14 .5) 

39.000 STOP 

40.000 END 















Dividing both sides of Eq. (5-24) by A and expanding the loga- 
rithmic term into a power series, we have 



ttA'gFO 


+ 


ttA^y 




\ 


i_ T 

ttA'gFO 


Y7T 


C-jA 




(5-32) 


Taking the limit on both sides of Eq. (5-32) as A-*0, and using the 
fact that 


tim C, = — — • 
A-K) 1 yT GF0 


(5-33) 


we have 


tim 

A-K) 



= tim 
A -^0 


T " 

?rA GFO 


Try 


C^A 





tim 

A-K) 


- — r + 

ttA'GFO 


4 V T GF0 , , 

TTY A ir 


Substituting Eq. (5-33) into the last equation we have 


J 


tim — L = tim +■ A^y A T B + 1 = o (5-' 

A+1 l A J A+0 GF0 Y GFO ■•• j 0 ls ' 

Dividing both sides of Eq. (5-30) by A and taking the limit as A 
approaches zero, we have 

fB,i ?f c i 1 

ti m I — j- 1 = tim 

A-0 ' A J A-*0 yA^ ^-A 2 


A "*° yA Vma/C, j 2 


(5-35) 


2 1 /2 

Expanding [1 - (A/Cj) ] ' into a power series, and using only the 
first two terms, Eq. (5-35) becomes 


? f 1 A 

Zcm = tun —y i- §- 

A-K) R i A-*0 yA 4 L 1 


Thus , 


Zajh\ o "* TT^pn 
A-+0 yC* 6 F0 


Um 1/N(A) = Um i — = -JL.-k - 

A-0 A*0 !i _ . V Y W 


(5-36) 


(5-37) 


As shown in Fig. 5-1, the gain-phase plot of -1/N(A) as A approaches 
zero is a point which lies on the -180 deg line with a magnitude of 

20 1°9]q^/yT GF q 3- 
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Asymptotic Behavior of -1/N(A) For Very large Values of A 
For very large values of A, the value of C-j becomes 


tim C 1 = tun 
A-*» 


A + 


YT 


6F0J 


(5-38) 


Then 


tun 

A-** 


l l _ 


tim 

A-+00 


4 2 

ttA ^GFO + Try 


c 7 a 


A 3 

3 


'1 U 1 


\ 

• *! 

t 

f 


= jj£! l" vrA T GF0, 


- -0 


Similarly we can show that 

tun -r— - +0 
A-*» a 


Thus, 


(5-39) 


(5-40) 


1/N(A)]» -j/0 °°q | -270 ( 5 . 41 ) 

As shown in Figure 5-1, the gain-phase plots of -1/N{A) approach 


> Ll?7P 0 as A •* » for all values of y and T™< 

GFO 


6. Computer Simulation of the Simplified LST System with the 
Analytical Torque Expressions 

A computer simulation of the LST system is presented here to 
corroborate the results of the describing function analysis of the 
last chapter. Since the describing function analysis has been carried 
out with the analytical torque expressions for the CMG frictional 
nonlinearity, the simulation model of the nonlinearity also has the 
same characteristics. This model of the nonlinearity is implemented by 
using the expressions for T QF in Eqs. (2-38) and (2-39) with initial 

conditions for 6 q and T^p being redefined each time a sign change in 

6 occurs. 

G 

The simplified LST system is represented by the block diagram of 

Figure 1-7. The linear transfer function which the nonlinear element NL 
sees is given by 


j g j v s4+j v k p s3+j v k i s2+k i hk 1 s+k i hk o 

Two sets of numerical values are considered as follows: 



System 1 

System 2 

J v 

10 s 

10 s 

J G . 

2.1 

3.7 

V 

216. 

280. 

K I 

9700. 

10000 

H 

600 

200 

K 1 

1371.02 

3000 

K o 

5758.35 

20000 


The f re quency- domain plots of G(s) for both systems are given 
in Figure 6-1 in db versus phase coordinates. Figure 6-1 also 
contains the -1/N curves of Figure 5-1 for y = 1 .38 x 10 5 , 

1.38 x 10 6 and 1.38 x 10 7 . 

With y = 1.38 x 10 7 , the -1/N curve intersects the G(s) curves 
of the two systems at two points each. Among these the stable points 
for sustained oscillations are the ones on the left at the higher 
frequencies. The approximate magnitudes and frequencies of the 
oscillations are 6 x 10 ” 6 rad and 4.4 rad/sec, respectively, for 
system 1, and 2 x ICf 5 rad and 5.6 rad/sec, respectively, for system 
2. The curves in Figure 6.1 also show that for y considerably smaller 

than 1.38 x 10 7 , both systems will exhibit a stable response, although 

\ 

for certain values of y system 2 will show sustained oscillations 
while system 1 is stable. 

For the computer simulation, the input to the LST system, x* 
is set to zero, along with all the initial states, except for the 
vehicle position e y . The initial value of e v is set at 5 x 1(T 5 
rad, which is chosen so that the input signal to the nonlinearity, 

0 g, would be large enough to cause the torque to saturate, while 
at the same time the limiting value of the input signal is not 
exceeded. 

The following quantities are plotted from the simulation runs: 

0 V = vehicle position (rad) 

= vehicle velocity (rad/sec) 

0 G -■ Gimbal position (rad) 













= Gimbal velocity (rad/sec) 

TqP s Torque output of the nonlinearity (ft-lb) 

Error = Error input commandt (rad/sec) to the CMG 

= X ~ ~ l<i<iJy 

Figures 6-2 and 6-3 show the plots of the above listed quantities 
for system 1 with y - 1.38 x 10 7 . It may be noted from the plot of 
in Figure 6-3 that the system has a sustained oscillation. This 
oscillation is not seen on the other plots because of the large initial 
transients. Figure 6-4 through 6-5 show the continuation of Figures 6-2 
and 6-3 with proper scales. Figures 6-6 and 6-7 show the response plots 
for system 2 with y =1.38 x 10 7 , and Figures 6-8 and 6-9 show the 

r 

continuation of these plots with proper scales. The frequencies and 
magnitude's of oscillations obtained with the two systems are quite 
close to the predicted values. The small discrepancy is attributed to 
the discretization of the nonlinearity implementation on the digital 
computer. ■■= 

Figures 6-10 and 6-11 show the response plots for system 1 with 
y- 0- 69 x 10 7 and their continuations are shown in Figures 6-12 and 
6-13, respectively. Figures 6-14 and 6-15 show the response plots 
for system 1 with y = 1.38 x 10 s . As predicted, the system is stable 
for both of the lower y values. 








O'O 



















7. 


Transfer Functions of the Sampled-Data LST System 


Since the actual LST system has sample-and-hold between the vehicle 
controller and the CMG controller, the system should be modelled as a 
sampled-data control system. Figure 7-1 shows the block diagram of the 
simplified LST system with sampled data. Since it is necessary to 
isolate the CMG nonlinearity from the linear dynamics for analytical 
purposes, a sample- and- ho Id is inserted in front of the nonlinearity 
as an approximation. 

Referring to Figure 7-1, the following equations are written 
using e*, e G * and 0y* as outputs. 


e* = (G ]X )* +N* 


rG ho G l G 6 G 7 


'0 


" *■ 

0 * 

" G ho G l G 2 G 3 G 6 G 7 

G 

A 0 i 


(7-1) 


0Q* = -N* 


G. G.l* 

--Q. 0*4. 

sA„ U G 


*0 J 


\ohW 


SA 


0 J 


e* 


(7-2) 


e v * = 


G* G«G~G^G-,1* 
ho 2 3 67 


e* - N* 


G. G.GVl* 
ho 6 7 


0 G* 


(7-3) 


where the symbol * denotes the z-transform operation, and N* represents 


the discrete describing function of the CMG frictional nonlinearity; 


A 0 ° 1 + G 3 G 6 (7-4) 

Equations (7-1) through (7-3) are portrayed by the sampled signal flow 

graph of Figure 7-2. Applying Mason's gain formula to this flow graph 
yields the determinant of the graph as 














Figure 7-2. The Sampled Signal Flow Graph for the Equivalent System 


81 


A - 1 + N* 


+ W* 


r 

“ G ho G 6l* [ G ho G 6 G 7 G ll 

* 

G ho G 2 G 3 G 6] 

*' 


sA o i L A) J 


sA o 

. 


’ G ho G l G 2 G 3 G 6 G 7 

* 

[ G hoV 

* 

4- 

r G ho G l G 2 G 3 G 6 G 7‘ 

■ A 0 J 


sA 0 J 

T 

*0 J 


(7-5) 


The last equation is put into the form of (with N(z) = N*) 
1 + N(z)G(z) = 0 


where 


G(z) = 


_ A 1 (z)-A 2 (z)A 3 (z)+A 4 (z)A 1 (z) 


A-j(z) = (1-z -1 


HA 4 (zT 

re, ) 


s 2 A, 


0 


-1 


A^z) = (1-z" ) £ 


A 


-1 


3 (z) « d-z" ) j 



-1 


A 4 (z) = <l-z~ ) ^ 


fG 1 G 2 G 3 G 6 G 7 1 


[ sA 0 J 


(7-6) 

(7-7) 

(7-8) 

(7-9) 

(7-10) 

(7-11) 


Substitution of the system transfer functions into the above expressions 


yields ? 

J„s>K s+K, 
A = -A ’.P 1 

0 ] s 2 

J Q s 


A,(i) = O-r -1 )}' 


1 


,J p s 3 +K s 2 +K t s, 

[ G P I J 


(7-12) 


(7-13) 


AgU) = (1-z" 1 ) 


H(K 0 +K n s) 


J v s 2 (J G s Z +K s+Kj) 


(7-14) 


A 3 (z) = (1-z' 1 ) j. 


K, 


-1 


A 4 (z) = (1-z *) J 


s 2 (0 G s2 +K p s+K I )j 

f KjHfl^s) | 


(7-15) 


J v s y (J G 7+K p s+K I ) 


(7-16) 


V 


The following system parameters are used for System 1 : 

H * 600 ft-lb-sec 
Jq * 2. 1 ft-lb-sec 2 
K q = 5.75835 x 10 3 
K 1 =1.37102 x 10 3 
K = 216 ft-lb/rad/sec 
Kj = 9700 ft- lb/ rad 
J v = 10 5 ft-lb-sec 2 

Taking the z-transforms of the functions Inside the brackets in 
Eqs. (7-13) through (7-16), we have the following results: 


(z) = A 13 + A 14 


z-1 . 

aT + A 15 


z-1 


z-e 


z-e 


aT 


(7-17) 


a = r 53.429993- j44. 42923 

a" = - 51.429993 + j44. 42923 

A ] 3 =1.0309376 x 10" 4 

A u = - 5.15469 x 10" 5 - J5. 9669183 x 10" 5 

A 15 = -5.15469 x 10' 5 + J5.9669T83 x 10" 5 


A 2 (z) " A 22 FT + ^3 + ^4 


Z-1 


t-e 


aT 


*25 


z-1 


z-e 


XT 


(7-18) 
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A 22 = 3.5616776 x 10" 3 

A 23 = 7.6870434 x 10' 4 

Ag 4 * -3.8435217 x 10" 4 - j4. 8499764 x 10 -4 

A 25 - 3.8435217 x 10' 4 + j4. 8499764 x 10" 4 

A 3 (z) = A 32 i^T + A 33 + A 34 ~af + A 35 ~aT 

Z“c Z “c 


a 3 2 1.0 

A 33 - -2.2268891 x 10‘ 2 

A 34 = 1.1134446 x 10' 2 + jl. 6350537 x 10" 3 

A 35 = 1.1134446 x 10' 2 - j 1.6350537 x 10" 3 


A 4 (z) = A 41 


T^(z+1 ) 
2 (z-1 ) 2 


A 42 z-1 + 


A 43 + 


y - z .-1 
44 z-e aT 


+ A 


A 41 = 34.54837 
A 42 = 7.4564409 
A 43 = -0.17352653 

A 44 =8.6763263 x 10' 2 + jl. 6520832 x 10' 2 
A 45 = 8.6763263 x 10' 2 - j 1.6520832 x 10' 2 

The following system parameters are used for System 2, The 
same expressions for A-j(z) , A 2 (z), A 3 (z) and A 4 (z) are preserved. 
H - 200 ft-lb-sec * 

Jg = 3.7 ft-lb-sec 2 
I . K 0 = 2 x 10 4 
, K 1 = 3 x 10 3 


(7-19) 


Z-1 

45 z-e 17 
(7-20) 


280 ft- lb/ rad/sec 
10 4 ft- lb/ rad 
TO 5 ft-Tb-sec 2 

The corresponding coefficients in Eqs. (7-17) through (7-20) 
are ' ' ; ■; 

a - -37.83783 - j35. 651077 
a = -37.83783 + j 35.651077 
A 13 * 9.9999976 x TO' 5 

A 14 = -4.9999973 x 10" 5 - j5. 3066848 x TO -5 

A, -= -4.9999973 x TO" 5 + j 5. 3066848 x TO" 5 

A 22 « 3.9999932 x TO* 3 

A 23 = 4.8799929 x TO' 4 

A^ = -2.43997 x 10' 4 - J3.T50655 x TO* 4 - 

Agg = -2.43997 x 10* 4 + j 3. T 50655 x TO' 4 < 

A 32 = 9.9999982 x TO" 1 

A 33 - -2.8 x 10' 2 

A 34 = 1. 4000002 x TO' 2 + j8. 3391 555 x TO' 4 

A 35 = T. 4000002 x TO' 2 - j8. 3391555 x 10' 4 

A 4 . = 39.999985 • v'' , . 

A 42 = 4.880001T 

A 43 = -T. 5144014 x 10' 1 

A 44 = 7.5720072 x TO' 2 + J1.T92342T x TO' 2 

A 45 = 7.5720072 x TO' 2 - jT. 1923421 x 10‘ 2 

It can be shown that if T approaches zero, the z-transfer function 
G(z) in Eq. (7-7) reverts to that of the continuous transfer function 
G(s) of Eq. (1-16). 


V 

K i = 

V 
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The G(z) plots for system 1 and 2 are shown in Figures 7-3- through 
7-6. Since these plots are to be used for predicting the existence of 
self-sustained oscillations in the sampled-data system, the frequency 
of oscillations is expressed as an integral multiple of the sampling 
period T. Thus, the frequency of oscillation is represented as 

2n (■ 

w c “ NT 

for N = 2, 3, . .. . 

The various plots of each system correspond to different values of 
N, with T as a parameter. Figures 7-3 and 7-5 are for lower values of 
T where as Figures 7-4 and 7-6 are for higher values of T. 

Notice that for large N and small T, the G(z) curves approach that 
of the continuous system transfer function G(s) . 
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8. The Discrete Describing Function of the CMG Frictional Nonlinearity 

In order to study the condition of self-sustained oscillations 
of the LST system with sampled data, it is necessary to evaluate 
the discrete describing function of the CMG frictional nonlinearity, 
N(z). 

The first step in the derivation of N(z) involves the inter- 
changing of the positions of the nonlinearity and the zero-order 
hold in Figure 8-la. This step is justified since the nonlinearity 
is amplitude dependent only, so that the signal of T^ is not affected 
by this interchange. Figure 8-lb illustrates the transposition 
between NL and zoh. 

The second step involves the assumption that 0^ is sinusoidal; 
that is , 

0q( t) = Acos(ojt + <f>) (8-1) 


where A, co, and denote the amplitude, the frequency in radians, 
and the phase in degrees of the sinusoid, respectively. 

The z-transform of 0g(t) is 

e G (z) « l Acos(^p + <j>)z _k (8-2) 

k-0 N 

or in closed form. 


0g(z) = 


p p 

A z[(z - cos|p-)cos4> - sin^sin^] 


jf 

- 2zcosj^- + 1 


(8-3) 



e G e G 

NL 

T * 
GF 

. zoh 

H ^ 

\ 


• 

L.~. — ■ — ^ 


(b) 

Figure 8-1 


An important consideration is that because of the periodic 
nature of the sampler, 6^(t), 0 G *(t), and ) are all periodic 

functions of period NT, where N is a positive integer >2. 
Therefore, to = 2tt/NT, and wT = 2 tt/N. 



The output of the nonlinearity in Figure 8-lb is denoted by 
T GF *(t), and its z-transform is Tgp(z). The discrete describing 
function (DDF) of the nonlinearity is defined as 


N(z) 


t gf (z) 

“Oil 


(8-4) 


It turns out that the discrete describing function (DDF) for 
N * 2 must be derived separately, and a general expression for N(z) 
can be obtained for all N > 3, 


The DDF for N = 2 

Let Tgp(kT) denote the value of T GF *(t) at t = kT. For N = 2 , 
the signal Tg F *(t) is a periodic function with a period of 2T. The 
z- transform of TQp*(t) is written 


T gf (z) ■ Tfip(O) (T + z~ 2 + z‘ 4 + ...) + T cc (T).(z" 1 ;+ z“ 3 ' + .. 


GF 


GF' 


T gf (0) z +■ Tgp(T)z 


2 . 
z - 1 


(8-5) 


For the CMG frictional nonlinearity, it has been established 
in chapter 2 that 


T GF (t) 


R 

t m 

GFQ 1 
R + 1 


■ g{1 -■ COS (lot + $)] 

+ ^{1 - cos(ut + <(>)] 


®G - 0 


(8-6) 
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R » a 


X ft N . T rT - T + 2 0 - C0S(Wt + ^) ] 

6F ' -«.«*•)] 


Let us introduce the following notation: 


T GF = T GF^'t*kf 0 G - 0 


T GF ^ “ T GF^ 


t=kT e G-° 


We have. 


T eF "(kt) = T 


R ar, 
R + 1" 2*- 1 


cos (^ + *)] 


GFG -A + gfi 
F+T 2r 


CO s(ir+ ^)] 


R .a. 


x + tm . T r ^ t + P - cos (ir + 

VKU - Jccri i “ ^fE 


e . 


GFO 1 a r , /2irk , m 
r“=“T + 2 il - cos(-|p + 4>)] 


For N = 2, ( 4> / ± tt/2) 


t G f(°) 


t gf ( t ) * 


V (0) 

0 _< $ < TT 

t gf + (°) 

TT < <f> < 2lT 

T gf + (T) 

0 < (p < 7]- 

yr 

i 

Ll. 
CD 
1 — 

tt < <£ < 2tt 


>0 (8-7) 

(8-8) 

(8-9) 

<0 ( 8 - 10 ) 
>0 ( 8 - 11 ) 

( 8 - 12 ) 
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Substituting N a 2 into Eq. (8-3), we have 


n Azcos<i> 

e G (z) - T+T 


(8-13) 


Using Eqs. (8-5) and (8-13), the DDF for H * 2 is determined. 


. yo). * y(T) 
1 } A(z - TJcosJ 


(8-14) 


Also, for N = 2, z = -1, the last equation becomes 


N(z) 


Tgp(O) " 


GF’ 


2Acos<j> 


(8-15) 


For stability analysis, we define 


F(z) a - 


2Acos^ 


nTzT ‘ t GF (tj - r GF (o) 


( 8 - 16 ) 


The DDF for H > 3 

In general, the z-transform of the output of the nonlinearity may 
be wri tten as 

^ • i. !'•**** 


N-l 

I Tgp(kT)z 

k=0 bh 


N-k 



(8-17) 


Using Eq. (8-2) for e G (z), the discrete describing 
N(z) is written 


T eF (z) 

>( '> • yrr 


N-l N Ir 

l T SF (kT)z N k 
k*0 w 


(z N - 1) 


k*0 




The denominator of N(z) may be simplified as follows: 


(z N - 1) l Acos(~^+ d>)z" k - A l z N " k cos(~p+ *) - A 
k=0 " k«0 N I 


N-l 

■ A l co s(^~ + <J>)z 
k=0 M 


2 irk , ^_N-k-l 


Thus, 


N(z) * 


N-7 .. . . 

I Tgc(kT)z^ - k " 
k=0 


aT7os(^ + *)z N - k -’ 

k-0 H 


(n > : 


As an alternative we may expand z N - 1 as 


.N 


N-l 


z - 1 = H (z - e 
k=0 


j2irk/N 


Then, using Eq. (8-3) for 0 G (z), N(z) is written 


function 


(8-18) 


I^ z k cos(—^ + d>) 
(8-19) 


( 8 - 20 ) 


( 8 - 21 ) 


( 8 - 22 ) 


N(z) 


lT r p(kT)z N " k [z 2 - 2zcosjjr + 1] 
k=0 N 

H (z - e^ 2ltk/,N )Az[(z - cosl^cosif) - sln^^i n<f>] 
k=0 " " 


For N « 3, N(z) Is simplified to 

. yo)z 2 ^ T gf (T)z + T gf (2 T) ' 
Nlz ' ' A(z - D[(z + Ojcos? - 0.866s1n<(>J 


(8-23) 


For N > 3, 


N(z) 


N-l 


N-k-1 


N-2 

a(z - 1) n (z 
, k«2 


e j27rk/N)[( z _ cos~)cos<^ - sin^si ncp] 


(8-24) 


where in general, 

T GF (kT) = T GF '(kT) 0<^+*<ir 
* Tg F + (kT) tt < ^jp- + <fi < 2 tt 


(8-25) 


where 0 < (2irk/N + ) < 2tt must be satisfied by appropriate conversion 

of the angle 2irk/N + 4 , 

For stability studies the critical regions of F(z) = -1/N(z) 
should be constructed for N = 2, 3, with <p varied from 0° to 
360°, and A from 0 to infinity. 


The following theorems on the properties of -1/N(z) are useful 
for sinplifylng the task of the construction of the critical regions. 


TktoAZjm &-1 

Voa any A.ntzgAaZ W, the. magnitude, and pkaAZ oh -7/W(z] a tpzat 
hon dvzAjy <j> - Zir/W AadlcunA. 

Proof: The negative inverse of the discrete describing function 

is written 


F(z) = - 


N (zl " 


■A[(z - cos~-)cos<fi - sin^sin1)](z N - 1) 


N-l 

1 ^Qp(^) z 

k=0 


N-k-1 


(z 2 - 2zcos§^ + 1) 


(8-26) 


Let 


Mz)F (z) 

F(2) = W 


(8-27) 


where 

F^ (z) - -A[(z - cos|^)cos<j) - si n^psi n<f>] 

F 2 (z) = z N - 1 

N-1 M t , 

F 3 (z) = l T G|r (kT)z^” k ' 

3 k=0 GF 

F 4 (z) = (z 2 - 2zcosj^- +1) 


(8-28) 

(8-29) 

(8-30) 

(8-31) 
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Let 


Then, 


and 


F(z) n = [F(z)] ( j i s ^ + 2lr /i\| 
F 1 (z) n = [F-j (z)]^ „ ^ + 27f y N 

f 3 (z)n = [F 3 (*)] t = <(, + 2 tt/N 



F-j (z) 

F 3^N 


F 1 (z) N 


(8-32) 


(8-33) 


(8-34) 


(8-35) 


Arg[F(z)] - Arg[F(z)^] = ArgtF^z)] - Arg[F 3 (z)] 

- Arg[F 1 (z) N l + Arg[F 3 (z) N ] (8-36) 

Also, 


F-j (z) 


Asinj^-' 


F 1 (2) N 


Asinj^ 


Let us express F 3 (z)j, as 



(8-37) 

(8-38) 


F 3^N = J 0 T GF^ kT) N zN 

where 

T GF^ kT ^N " ^ T GF^ kT ^<j) = <t> + 2-rr/N 
It can be shown that for any integral N, 

F 3^N " zF 3( z > = " = V GF(kT)zN ' k 
Then, Eq. (8-35) becomes 


F 

1 

J 

F 3^N 

i 1 ZF 3 (Z) 

F 

[z) 

N 



\ W 


The argument of F^(z) is 

Arg[F^(z)] = <)> - ir/2 


Then 


Arg[F-,(z) N ] = <(> - n/2 + 2tt/N 


Thus, Eq. (8-36) becomes 


(8-39) 


(8-40) 


(8-41) 


(8-42) 


(8-43) 


(8-44) 



TOO 


Arg[F(z)] - Arg[F(z) N ] » - jp- - 



2tt 

" ' N ' 


Arg[F 3 (z)] + Arg[F 3 (z) N ] 
Arg[F g (z)] + Arg[zF 3 (z)] 
Arg[F 3 (z)] + jp+ Arg[F 3 (z)] 


Q.E.D. 


*• 0 


As an illustrative example of Theorem 8-1, let us consider the 
case of N - 4. 

Let | = Y T gfo A = K 

Then 


t gf ( kT ) s 
t gf = 



For 0 < <J> < tt. 


F 3 (z) = T gf '(0)z 3 + T gf '(T)z 2 + T gf + (2T)z + T gf + ( 3T) 
F 3^N = T GF ^°^N z3 + t gf + t gf < 2t V + t gf ( 3T >n 


It is easy to see that 


(8-45) 

(8-46) 

(8-47) 

(8-48) 


(8-49) 
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T GF = T GF ^ 

T GF < T >N = T GF ( 2T ) 
t gf < 2t > N = T GF ( 3T ) 

T Gf'< 3T V = V (0) 

which proves that F 3 (z)^ - zF 3 (z). 

Similar results are obtained for tt < <J> < 2tt. 

Tho.oh.zxn 8-2 

Foh odd N (A/ > 3), tkz magnitude, and phabz oj 5 - / /W t z ) hz.pz.at 
&oh zvzhy <p * tt/A/. 

Proof: Let F(z)^, F-j(z)^, F 3 (z)^ now be defined as F(z), 

F-j(z), and F 3 (z) with ^replaced by <t> + ir/H, respectively. 
Then, 

Arg[F 1 (z) N ] = <j> - 2- 

\ 

Arg[F(z)] - Arg[F(z)^j = - j|- + Arg[F 3 (z) H ] - Arg[F 3 (z)] 

Using the same notation as in Theorem 8-1, it can be shown that 
for odd N ^ 3, 


(8-50) 


(8-51) 

(8-52) 


102 


F 3 (z) n - -z- (N - 1)/2 F 3 (z) (8-53) 

Thus, 

Arg[F 3 (z) N ] = tt - + Ar g[ F 3( z )] 

= ^ + Arg[F 3 (z)] 

Again, 

i- 

■ 1 
and 

Arg[F(z)] - Arg[F(z) N ] = - f + Arg[F 3 (z) N ] - Arg[F 3 (z)] 

=0 (8-56) 

Q.E.D. 

As an illustrative example of Theorem 8-2, consider the case N * 3. 

For 0 < <j> < it > 

F 3 (z) = T gf '(0)z 2 + T gf ‘(T)z + T gf + (2T) (8-57) 

>V Z >N = + T GF + < T ) N Z + T GF + ( 2T) N (8-58) 

It can readily be shown that 



(8-54) 


(8-55) 
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V< 0 >n = V( 2T > 

t gf +(t) n = ' t gf'(°> (8 - 59) 

t gf + ^n = -V^ 

To carry out one of the identities above, 

R - l)(cos<j> + 1) T 
R -H)(cosi + i)'gfo 

(8-60) 
(8-61) 

Using the relation, R 2 + R/K - 1 = 0, we have 

T GF' (T) N * 'V (0) (S-62) 

Thus, 

F 3 (z) n » -z' 1 F 3 (z) (8-63) 

The significance of the last two theorems is that the critical 
regions of -1/N(z) need be computed only for 0° < <f> < tt/N for odd 
N, and 0° < $ _< 2ir/N for even N. 


GF 


’ T GF0~1 ‘ 


- K(R - l)[cos($ + tt) - 1] _ R + K( 

:n<TR— ijtcosU + ir) -'ll rri<t 


GF 


m\ _ -r R + K(R + l)(cos$ - 1) 
(0) ' T GF0~ 1 — - K(R +’ l)(C0Sij> - 1) 


The.otie.rn S-3. 

A&ympotic Behavior o£ -7/N(zJ at A appsioacheA infinity. 

(a) I - 1 /N (z ) | = « (8-64) 

A-*» 

(b) For even N >_ 4, 0 < 4> < 2 tt/N. 

tim Arg[-1/N(z)] = (i - [Sr) 11 + c|> . (8-65) 

For odd N > 3, 0 < <j> £ 7r/N. 

^ Arg[-1/N(z)] = (1 - i)f+ <f> (8-66) 

A-«= N ^ 

(c) For N " 2 0 <^£tt 

tim Arg[-1/N(z)] - 0° 0 < 4> < tt/2 

A-h» 

2 tt tt/ 2 < <j> < tt (8-67) 

(d) For N * 3 0 < £ tt/3 

Arg[-1/N(z)] = (8-68) 

A-*« 0 

Proof: 

We can easily show that 

* 

tim Tgp (kT) = ~Tqpq 9q <. 0 (8-69) 

k+co 


tim T 
A-*°° 


+ 

GF 


(kT) 





(8-70) 


The magnitude of -1/N(z) is directly proportional to A as A 
approaches infinity; thus (a) is proved. 

For N ~ 2 0 < <f> < tt/2 


Thus , 


F <-> ■ -wi ■ i. F ( T f“L(ff 


tim Arg[F(z)] - ton Arg ■■ 

Ar*» VgF UJ 'GF' 0 ', 


- Arg |^osi a Q o 0 < < tt/2 

GFO 

= IT tt/2 < $ <. IT 

This proves item (c). 

For N = 3 0 < <J> 5 tt/3 


F( z ) = -A[(z + 0-5)cos<j> - 0.866sin(ji](z - 1) 
T gf (0)z 2 + T gf (T)z + T gf (2T) 


tun Arg[F(z)] » - - Um Arg[T rF (0)z 2 + T rl -(T)z + T„ C (2T)] + <1 

A-x» ■ A-*» bh 6h 


j- - Arg[-z 2 


- z + 1] + ^ 


5ir 


+ <l> 


This proves item (d). 


(8-71) 


(8-72) 


(8-73) 


(8-74) 


For N >_ 4 and even 0 < <j> < 2 tt/N 

Using Eqs. (8-27) through (8-31), we have 


tun Arg[F(z)] * tun Arg[Fj(z)] + Arg[F 2 (z)] - Arg[F 3 (z)] - Arg[F 4 (z)] 
A~*» Pr*co V 


f - ^ [Arg[F 3 (z)]J + Arg[F 2 (z)] - Arg[F 4 (z)] 


(8-75) 


N-2 


J2irk/N 


Arg[F 2 (z)] - Arg[F 4 (z)] ■ Arg[(z - 1) IT (z - e JC,, * /M )] 

* k*2 


N-2 


- Arg[e J ^(e^ N - e~^ N ) H - e J 27r /N(k-l)^ (8-76) 

k*2 


Arg[F 2 (z)] - Arg[F 4 (z)} * Arg[e J7r ^(e^ N _ e -“J 1 >‘/Njj e j27r/N.jN-3 ! n ^ _ e^ 27rk//N )] 

k=l 


N-3 J 


- J + f +' (N - -3)(|3L) + Arg H y* k/ V~ jWN - e + ^ k/N ) 


k*l 


N-3, 


+ (N - 3)({*) ■+ (N - 4)(-|) + 'J,W 


k*l 


(8-77) 


N-l 

Urn Arg[F-(z)] * tun Arg J T rP (kT)z 
A-*° J A^» k=0 


N-k-1 
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w$ H-k-i N ; 7 N-k-i 

Arg I - T sfo z N k 1 + Arg l T GF / k 1 

k=0 wo k=N/2 


N/2-1 

I 

k=0 


m'.l- W 1 ' 1 '' 1 


if-')* 


ir 7T 
N “ 2 " N 


Thus , 


tim Arg[F(z)] 
A** 50 


♦ - f + f + (N - 4) (- j) + (N - 3)(f) + Vj* 


_ /o N. 4 . k\ . . 

3 " 2 " N k =f N * 


s ^ + $ 


For N > 3 and odd 0 < <p <_ tt/N 

i 

1 For this case, 


(N- 1 V 2 

tan Arg[F~(z)] = Arg J 
A-*» k=*0 


- T„ m z N ' k * 1 + ArgV 


GFO 


k$<+T)/Z 


t gfo z 


N-k-1 


. (N - 3)it 
2N 


(8-78) 


7T , 7T 

2 + N 


(8-79) 


(8-80) 


Thus, 


tcm Arg[F(z)] = <l>' ' jr + ft + (N - 4) (- ft) + (N - 3)(^L) + F" 

and (b) is proved. Q.E.D. 

Tkconjm S-4. 

Asymptotic Bckavtor. Ofj -7/W(z) as A approaches zvio. 
til n F(z) » W 

*° ^ t gfo 


for all <p and all N. 

Proof: From Eq. (8-20), for N 3, 


F(z) 


-mh = 


‘aVco»(2Ji+ *)z M ' k - 1 

k=0 N 

"i'vn)^' 

k=0 bh 


tan F(z) ■ 
A+0 


k=0 N 


tan \ 
A-*0 k=0 


N=1 T GF (kT) 


N-k-1 


(N - 3)ir 
2N 


(8-81) 


(8-82) 


(8-83) 


(8-84) 


N-l 


- I cos( 
k=0 


2kir 

N 


+ <J.)z 


' N-l 
l Um 
k=0 A+0 


T^kiy 


N-k-1 


N-k-1 


T GF (kT) 

Therefore, the problem is that of finding tim — - 

M A 

First, let Tg F (kT) = T GF "(kT), Then 


tim 

A-*0 


t sf (kT) T Arrrry ~ yT gfo + yT gfo cos( it + ^ 

A ™ GF0 vW • vW c ° s (¥ + *> + rh 


s T GF0 ^2 Al RTT ) + ^gfo 2[cos( T + ♦> - < 8 - 85 ) 


where the fact that [1/(R + 1)] ~ 1 has been used. 

A-*0 


^”o A(R+ n 


tim ~ = tim i 
A->0 A 


2yT GF0 A 


/4 y 2 T 2 a 2 + 1 

' GFO A 1 


4y T 


GFO 


2 A 2 


(8-86) 


Or, 


R 


^KFTtr 


- tim 
A+0 


2yT GF0 A 2 yT GF0 A 


AT 


4y 2 T 2 A 2 
n i GF0 a 


(8-87) 



I 

no 


r — — 2 FT 

Expanding*/! + 4y T rFn A into a power series, and using only the 


GFO 

first two terms, we have. 


w aif+tt 


lb 71 i 

A-H) * 


4y^T 2 A 2> * 
i GFQ a 


4yT GF0 A 


-yT 


GFO 


(8-88) 


Thus, 


lim 

A*0 


T GF 


y "^GF0 + y "^GF0 I- cos ^ w + 40 " ^ 


N 


= yy GF0 cos ^ n + 4>) 


(8-89) 


Similarly, it can be shown that 


tim 

A-*0 


T GF + < kT > _ V< kT > 


= lim 
A-*0 


A 


YTqpq cos( ^ + <f>) 


(8-90) 


Now, 


- I cos(^-+ *)z 


N-k-1 


lun\ F(z) 
A+0 


k=0 


2N T , eosfM+ *i 7 N-k-l 
GFO k | 0 ( N * >z 


YT 


Yt 


GFO 


(8-91) 


for all $ and all N >3. 




9. Discrete Describing Function Plots of the CMG Frictional Nonlinearity 

The discrete describing function, N(z), for the CMG frictional 
nonlinearity is derived in Chapter 8 for integral values of N ^2. 

The plots of -1/N(z) together with that of G(z) in the frequency -domain 
allow the study of the condition of self-sustained oscillations of the 
LST system with sampled data. 

For N - 2, the expression for F(z) = -1/N{z) is given by Eq. (8-16). 
Figure 9-1 shows the F(z) plot for N - 2 in the gain-phase 1 coordinates 
with 0 < A < * and various values of The value of y is 1.38x10^. Note 
that the plot stays on the -180° and -'360° axes. 

In general, N(z) for N > 3 is given by Eq.. (8-20). Figure 9-2 shows 
the gain-phase plot for F(z) when N - 3. The curves for several values of 
$ are plotted to illustrate the effect of varying the phase of the input 
signal to the nonlinearity. It should be noted that the values of F(z) 
repeat every 60 degrees starting from <f> = 0°. As the magnitude of the 
input signal, A, approaches infinity, the bounds of F(z) are at -240° 
and -300°. Figures 9-3 and 9-4 illustrate the F(z) plots for N s 4 and 
N - 5, respectively. For N - 4, the F(z) plot extends from -315° to 
-225°, and for N = 5, the span is from -288° to -252°. 

For stability analysis, it is sufficient to consider only the 
bounds of the F(z) plot for a fixed N. Self-sustained oscillations 
characterized by N may occur if G(z) intersects with any part of the 
F(z) plot. The region bounded by all the F(z) curves for a given N is 
defined as the “ critical region ". In Figure 9-5 the critical region for 
N - 6 is shown, without showing all the curves of F(z) for various 
values of <p, Similarly, Figures 9-6 through 9-8 show the critical regions 
for N - 8, 20, and 50, respectively. The general shape of the critical 
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regions for other values of N is easily visualized. Furthermore, Theorems 
8-3 and 8-4 on the asymptotic behavior of F(z) as A *>'» and A o are useful 
in generating the critical regions. It is interesting to note that as N 
approaches infinity, F(z) approaches -1/N of the continuous-data non- 
linearity, as shown in Figure 9-8. 




Figure 9-1. Discrete describing function plot of CM6 frictional 


nonlinearity, y = 1.38 x 10 5 , N =2. 
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10 . 


Computer Simulation of the Simplified Samp led- Data LST System 
with the Analytical CMG Frictional Torque Expressions 

A computer simulation of the Sampled- Data LST -System is presented 
here to corroborate the results of the discrete describing function ' 
analysis of the last two chapters. Since the analysis has been carried 
put with the analytical torque expressions for the CMG frictional 
nonlinearity, the simulation model of the nonlinearity also has the 
same characteristics. 

The simplified sampled-data LST system is represented by the block 

diagram of Figure 7-1. Although for analytical convenience the 

, % 

discrete describing function analysis has been carried out with two 
samplers present in the system, the actual system has only one sampler. 
This is at the input to the CMG gimbal drive and there is no sampler 
in the nonlinearity loop. Consequently, the simulations have been 
performed with the two-sampler as well as the one-sampler system models. 
The numerical values used in the simulation model correspond to those 
of System 1 in Chapter 6. 

The G(z) plots of Figures 7-3 and 7-4 and the-l/N(z) plots of 

i E 

Figures 9-1 through 9-9 show that for y = 1.38x10 self-sustained 
oscillations will not exist if T is less than 0.25 seconds approximately 
For larger values of y the plots in Figures 9-1 through 9-9 shift 
downward appropriately. Thus, with y = 1.38xl0 7 the lowest point in 
the -1/N(z) curves becomes approximately -102 db, yielding a system in 
which self-sustained oscillations will always exist for any sampling 
period. The periods of these oscillations depend on the sampling 
period used. 
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For the computer simulation, the input to the LST system, x> is 
set to zero, along with all the initial states, except for the vehicle 
position 9y. The following quantities are plotted from the simulation 
runs 

9y « vehicle position (radians) 

= vehicle velocity (radians/second) 

- Gimbal position (radians) 

Wq = Gimbal velocity (radians/second) 

Tqf “ Torque output of the nonlinearity (ft-lb) 

Error = Error input command to the CMG (radians/second) 

x '^0 6 V : 

Figures 10-1 through 10-8 show the simulation results with y = 1.38x10^. 
Two different sampling periods, T * 0.005 sec and T = 0.1 sec are 
considered with both the one-sampler and the two-sampler system models. 

Table 10-1 provides the details of each of these simulation runs. Note 
that, as predicted, both the system models are stable with the sampling 
periods considered and no self-sustained oscillations exist. 

Table 10-1. Computer Simulations with y = 1.38x10 s . 

Figure No. No. of Samplers in Sampling Period . Initial Value 
System Model ft sec) ... of Gy (rad) 


10-1, 10-2 

2 

0.005 

l.xio* 

10-3, 10-4 

1 

0.005 

l.xlO' 

10-5, 10-6 

2 

0.1 

l.xlO' 

10-7, 10-8 

1 

0.1 

l.xlO- 


Figures 10-9 through 10-30 show the simulation-results with 
Y - 1.38x10^. Several sampling periods from T = 0.005 sec up to 
T - 0.25 sec are used with both the one- and 2-sampler system models 
The details of the simulation parameters and the oscillation periods 
are shown in Table 10-2. Note that the one-sampler and two-sampler 
system models alweys oscillate with periods close to those which are 
predicted by the -1/N(z) and G(z) plots. Although the two-sampler 
model is in closer agreement with the theoretical predictions, these 
results do justify the approximation of introducing the sampler. 

Note also that for some sampling periods (eg. 0.02 isec) although the 
one-sampler model and the two-sampler model oscillate with periods 
which differ considerably (0.04 sec and 1.76 sec, respectively) both 
periods are in fact predicted by the theory. 

From the standpoint of self-sustained oscillations, if y is at its 
nominal value of 1.38*10^ any sampling period less than 0.25 seconds, 
approximately should yield a stable system. However, other practical 
considerations and stability considerations due to non-zero Inputs 
would limit the sampling period to a much lower value. 
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\ . 

Table 10-2. Computer Simulations with y = 1.38x10 7 

i 

i . 


Figure No. 

No. of Samplers 
in System Model 

Sampling 

Period 

T(sec) 

Initial Value 
of 0y(rad) 

Osci llation 
Period (sec) 
(NT) 

.10-9, 

10-10 

2 

0.005 

l.xlO' 8 

1.58 

(316T) 

10-11, 

10-12 

1 

0.005 

l.xlO' 8 

1.88 

(376T ) 

10-13, 

10-14 

2 

0.02 

l.xlO' 7 

0.04 

(2T ) 

10-15, 

10-16 

2 

0.02 

Continuation 
of 10-13, 
10-14 

0.04 

(2T) 

10-17, 

10-18 

} 

0.02 

5.xl0~ s : 

1.76 

<88T) 

10-19, 

10-20 t 

2 

0.05 

l.xlO' 8 

0.1 

(2T) 

10-21, 

10-22 

1 

0.05 

l.xlO' 8 

1.25 

(25T) 

10-23, 

10-24 

2 

v 

0.1 

l.xlO- 8 

1.0 

(10T ) 

10-25, 

10-26 

1 

0.1 

l.xlO' 8 

1,2 

(12T ) 

10-27, 

10-28 

2 

0.25 

l.xlO' 7 

0.5 

(2T) 

10-29, 

10-30 

1 

0.25 

l.xlO- 7 

0.75 

(3T } 
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